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We propose an explanation of the superconducting transitions discovered in the heavy fermion 
superlattices by Mizukami et al. (Nature Physics 7, 849 (2011)) in terms of Berezinskii-Kosterlitz- 
Thouless transition. We observe that the effective mass mismatch between the heavy fermion 
superconductor and the normal metal regions provides an effective barrier that enables quasi 2D 
superconductivity in such systems. We show that the resistivity data, both with and without 
magnetic field, are consistent with BKT transition. Furthermore, we study the influence of a nearby 
magnetic quantum critical point on the vortex system, and find that the vortex core energy can be 
significantly reduced due to magnetic fluctuations. Further reduction of the gap with decreasing 
number of layers is understood as a result of pair breaking effect of Yb ions at the interface. 



Thin film growth technology recently has advanced to 
the point that artificial two-dimensional structures can 
be fabricated with atomic-layer precision. This has en- 
abled the exploration of novel aspects of emergent phe- 
nomena in low dimensional systems with unprecedented 
control. Using the molecular beam epitaxy (MBE) tech- 
nique, Mizukami et al. have grown CeCoIn5/YbCoIn5 
superlattices, where superconductivity was found to oc- 
cur in the two-dimensional Kondo lattice pQ. The com- 
bination of f-electron physics, low dimensionality and in- 
terface effects provides a rare opportunity to study new 
states in strongly correlated electron systems, e.g. uncon- 
ventional superconductivity, dimensionally-tuned quan- 
tum criticality 2 J , interplay of magnetism and supercon- 
ductivity, Fulde-Ferrell-Larkin-Ovchinnikov phases, and 
to induce symmetry breaking not available in the bulk 
like locally broken inversion symmetry [3] . 

Here, we investigate the mechanism for the onset of su- 
perconductivity in such heavy fermion superlattices. We 
propose an explanation of the experimental results of [T] 
within the framework of Bcrezinskii-Kosterlitz-Thouless 
(BKT) transition, and further study the interplay of 
Kondo lattice physics and BKT mechanism. While well 
established for superfluid films, BKT transition is less 
convincing for superconductors (See ?4 and references 
therein) . Though implications have been found in numer- 
ous thin superconducting films 4 8 , highly anisotropic 
cuprates [5HTT] , oxide interfaces [P2"HT4"] , the results have 
remained inconclusive (see e.g. |15ll6j ). It is therefore 
desirable to have a well-controlled, readily-tunable sys- 
tem to investigate the BKT physics. The epitaxially 
grown heavy fermion superlattices may serve such a role. 

Quasi 2- dimensional superconductivity: First, we dis- 
cuss why BKT theory is applicable to heavy fermion 
superlattices. In the CeCoIn 5 /YbCoIn 5 superlat- 
tice, one has a layered structure of alternating heavy 
fermion superconductor (CeCoIns) and conventional 
metal (YbCoIn 5 ), typically 3.5 nm thick. Prox- 
imity effect is expected to happen in such normal 
metal/superconductor (N/S) junctions. For conventional 
superconductors, the thickness of the leakage region is on 
the order of the thermal length Hvn /2Trk bT ', where vn is 



the Fermi velocity in the N region (see e.g. [T7] ). At low 
temperatures, this thickness is typically of order lOOrtm, 
which is much larger than the separation of CeCoIn5 lay- 
ers. One may thus expect a strong coupling between the 
superconducting CeCoIn 5 layers and the system would 
behave as three dimensional superconductor. However, 
as we will argue below, the large mismatch of Fermi ve- 
locities across the interface changes the story completely 
and enables quasi 2D superconductivity in CeCoIns thin 
layers. 

In normal metal/heavy fermion superconductor prox- 
imity effect studies, it was realized that the large mis- 
match of effective mass at the interface leads to huge 
suppression of transmission of electron probability cur- 
rents |18j . The ratio r<r of the transmitted probability 
current and the incident current is determined by the ra- 
tio of the effective masses, tt — imi/mh, for to/, 3> mi 
|18j . The effective mass of CeCoIns is of order 100m e . 
For the more conventional metal YbCoIn5, we take its 
effect mass to be of order m e . The transmission is thus 
on the order of one percent. 

This result is intimately related to that of Blonder, 
Tinkham and Klapwijk |19l20j . where it was shown that 
the mismatch of Fermi velocities between the N and 
S regions increases the barrier height between the two, 
with the effective barrier parameter Z modified to Z = 
(Zq + (1 — r) 2 /4r) 1 / 2 where r = vs/v^ is the ratio of two 
Fermi velocities. This gives essentially the same result as 
Ref. [18]. This suppression factor significantly degrades 
the proximity coupling to the point where 4 nm normal 
layer renders heavy fermion films essentially uncoupled. 
A direct consequence of the reduced proximity effect is 
an enhanced c axis resistivity, which can be measured 
directly in experiment. 

More extensive numerical studies of proximity effect in 
N/S junctions have been carried out recently [21] . where 
it was shown that proximity effect is substantially sup- 
pressed with moderate mismatch of Fermi energies. An- 
other source of suppression of the proximity effect is the 
pair breaking effects of Yb ions at the interface (see sup- 
plementary material). It is also expected that a weak 
magnetic field can destroy the proximity-induced super- 
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FIG. 1: Gap and T c as function of number of CeCoIns layers 
n (data from Mizukami et al. PQ). For n > 5 (shaded region), 
gap retains the bulk value, while T c decreases with decreasing 
number of layers. 

conductivity in YbCoIns layers |l|22j . 

Suppression of the proximity effect in the 
CeCoIns/YbCoIns superlattice and the fact that 
the thickness of the CeCoIns layers is on the order of 
the perpendicular coherence length £j_ ~ 20A pQ, lead 
to the conclusion that superconductivity in such systems 
is essentially two dimensional, and one expects BKT 
physics to be relevant in such systems. 

BKT transition: The basic experimental fact of 
Mizukami et.al [T] is that when the number of CeCoIns 
layers n > 5, the upper critical field H C 2, both paral- 
lel and perpendicular to the ab-plane, retains the bulk 
value, while the transition temperature T c decreases 
with decreasing n (see Fig[T]). H C 2 in such systems is 
Pauli-limited in both parallel and perpendicular direc- 
tions }1I23[ and is thus a direct measure of the super- 
conducting gap, with H^ 2 auh — VSA/g^B, where g is 
the gyromagnetic factor and ps is the Bohr magne- 
ton. This means that gap retains the bulk value for 
n > 5. The behavior of gap and T c for different num- 
ber of CeCoIn5 layers is shown in Fig. [I] Our proposal 
is that such behavior is due to the effect of phase fluc- 
tuations, which for the quasi-two-dimensional supercon- 
ductors considered here is controlled by the Berezinskii- 
Kosterlitz-Thouless physics [24125] . 

For two dimensional systems with continuous Abelian 
symmetry, despite the lack of broken symmetry due 
to strong fluctuations, there exists a finite temperature 
phase transition mediated by topological defects, e.g. 
vortices for superconductors [24 25 . Below the tran- 
sition temperature Tbkt, vortices and antivortices are 
bound into pairs, and the resistance vanishes. Above 
Tbkt, vortex-antivortex pairs unbind, and the prolifera- 
tion of free vortices destroys superconductivity. For such 
systems, one thus has T c = Tbkt- 

For layered superconductors, one also needs to include 
interlayer couplings. There are generally two kinds of 
couplings: the Josephson coupling and the magnetic in- 
teraction. Since the separation of the different CeCoIns 
layers is larger than the perpendicular coherence length, 
the interlayer Josephson coupling is weak, and can be ig- 
nored. The long range magnetic interaction couples vor- 
tices in different planes, and aligns vortices of the same 
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FIG. 2: Resistivity as function of temperature for n — 4, 5, 7, 9 
(data from Mizukami et al. [T]). The BKT transition tem- 
perature is determined from the intersection with the T-axis 
to be Tbkt = 1.202, 1.344, 1.582, 1.712if respectively. 
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FIG. 3: The BKT transition temperature Tbkt as function 
of the number of CeCoIns layers. The dashed line is T c q = 
2.3K. The solid line is a fit to the theoretical result, with 
I=p |ta_ F= 0.444(n- 0.317). 



sign into stacks. Since the interlayer coupling is still loga- 
rithmic as in two dimensional superconductors, the phase 
transition is expected to remain in the same universality 
class as BKT transition ,26] . This has been confirmed by 
detailed renormalization group studies [27H30] (see also 
[31] ). It has also been shown in Ref. [3D] that T c is only 
slightly modified.^ While such small modification may 
be detected by future high precision measurements, as 
first approximation we will ignore it in the following and 
concentrate on the single-layer problem. 

In the following, we are going to check whether the ex- 
perimental findings of Mizukami et al. [1] are consistent 
with BKT transition, i) First, we will examine whether 
resistivity has the right temperature dependence, ii) 
Then we extract from the resistivity data the transition 
temperature Tbkt- iii) Finally, we will check whether 
Tbkt has the right dependence on the number of layers. 
We find that the observations in [T] are consistent with 
BKT transition. 

Near T B kt, resistivity behaves as p{T) = 
/ 9oe- b(T - TBKT) " 1/2 [3J, which gives (din p(T)/dTy 2/3 = 
(2/bf /3 (T - Tbkt)- We plot in Fig. ^ the temperature 
dependence of (d In p{T) / dT)" 2 ^ for the four different 
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cases with number of CeCoIri5 layers n = 4,5, 7, 9, where 
one can see that (d\np(T)/dT)~ 2 / 3 is indeed linear in T, 
and Tbkt can be extracted from the intersection points. 
We also notice that resistivity does not fall to zero at 
Tbkt- It retains a small nonzero value in a temperature 
region below Tbkt- This is generically observed for a 
BKT transition, and is attributed to the temperature 
difference between the formation of single vortices and 
the subsequent vortex condensation (see e.g. and 
references therein). Inhomogeneity and finite size effects 
also broaden the BKT transition, giving rise to the 
resistivity tail below Tbkt |16j . 

Now, we proceed to study the thickness dependence of 
the BKT transition temperature. Tbkt can be written 
as [^51321341351 
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with the dielectric constant e c 
is the renormalized carrier density. The unrenormal- 
ized 2d carrier density n 2D — n 3D d is determined by 
the 3d carrier density n 3 s D (T) = n^(0)A 2 (0)/A§(T), 
and the film thickness d. The bulk penetration depth 
Af,(T) has a temperature dependence of the form \%{T) — 

A b (0) [1 - (T/T c0 ) Q ] _1/2 , with bulk mean field transition 
temperature T c q. In the usual two-fluid picture, the ex- 
ponent a = 4. For cuprates and CeCoIns, it has been 
found that a = 2 36 37J. Thus we have 



Tbkt 



1 - (T B kt/T c0 ) 5 
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Noting that d = nx — do = (n — no)x, with n the number 
of CeCoIns layers, x the thickness of each layer and do 
the thickness of the dead layers on top and bottom, the 
above result can be written as 
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We plot in Fig. [3] Tbkt as function of the number of 
CeCoIn 5 layers. The experimental results are in good 
agreement with the theoretical prediction determined 
from Eq. § Taking A 6 (0) = 358nm [37], x = £ c /4 = 
2.1nm/4, we get the fitting parameter e c ~ 90. With 
A~ 2 = XT /e c , our prediction is that the penetration 
depth of the superlattice is enhanced by about one order 
of magnitude from the bulk value. Furthermore, another 
important prediction from BKT transition that can be 
checked is that the penetration depth of the superlattice 
A satisfies the universal relation EH 
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right below the transition temperature, where $o = 
hc/2e is the flux quantum. 

Antiferromagnetic vortex core: We extract from the 
experiment [TJ a large dielectric constant e c , which in- 
dicates a large fugacity, or a small vortex core energy 





FIG. 4: Change of vortex core energy as function of distance 
to the QCP and the resulting a — T phase diagram as deduced 
from simulation [42] and theoretical results [43144] (see also 
[45146] h Phase A is gas of free vortices, B a gas of bound 
vortex-antivortex pairs, C a crystal of vortices and antivor- 
tices, D a hexatic phase of vortices and antivortices. 



[25134] (see supplementary material for a more detailed 
analysis). Here, we try to understand where such a large 
renormalization may come from. We find that at the vor- 
tex core, where the superconducting gap is suppressed, 
magnetic ordering can occur locally (see e.g. [38] h which 
reduces the vortex core energy. 

A salient feature of the heavy-fermion superconduc- 
tor CeCoIns is the proximity to an antiferromagnetic 
quantum critical point (QCP). Therefore, one may ex- 
pect that fluctuating magnetic order may influence the 
vortex dynamics in the heavy fermion superlattices. Sup- 
pression of the superconductivity in the core can induce 
the antiferromagnetic state in the cores as opposed to a 
simple metal in conventional superconductors. To model 
this effect, we consider magnetic moment that couples 
to the vortex via a Zeeman term g[i B HlS z , where H* 
is the magnetic field generated by vortices. H% is a 
superpostion of the magnetic fields generated by vor- 
tices at different locations, H*{r) — Y^i n iHo{ r ~ Ri), 
with rii the vorticity. Ho(v) can be obtained from its 
Fourier transform Ho(k) = $0/(1 + A 2 k 2 ), with re- 
sult H (r) ~ ($ /A 2 )ifo(r/A), where K is the mod- 
ified Bessel function of the second kind. For r <C A, 
K (r/X) ~ lnr. 

Zeeman coupling induces a precession of the magnetic 
moment perpendicular to the magnetic field, which can 
be captured by modifying the kinetic energy density to 
(d T cf) + igfj, B H x cj)) 2 , where </> is the sublattice magneti- 
zation density [55H4T] . For H in the z-direction, one can 
define $ = {<p x + i(j) y )/\/2. Consider the static limit, its 
free energy density reads 

T* = |V$| 2 + (a - g 2 v 2 B H 2 (r))\<P\ 2 + 7 |$| 4 



(5) 



Near the vortex core, H ~ In |r — r,| can be very large. 
Close to the QCP, a is small. When a = a—g 2 (i B H 2 < 0, 
the vortex core becomes antiferromagnetic, and quali- 
tatively |<I>| 2 = — d/27 and the potential energy V$ = 
— d 2 /47 < 0. Thus the vortex core energy is significantly 
reduced due to magnetic fluctuations. 

More precisely, we consider the equation of motion 

(_V 2 + a - g 2 ^ 2 B H 2 (v) + 2 7 |$(r)| 2 )$(r) = 0, (6) 
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where a vortex of unit vorticity is placed at r = 0. Far 
away from the vortex core, i.e. r > A, Hq decays expo- 
nentially, and $ = is the lowest energy solution. Near 
the vortex core, we can ignore a and $(r) ~ ln(r/A) is 
the lowest energy solution. The change of vortex core en- 
ergy is SE C = f d 2 r.F[$(r)] ~ - 5 4 /4*o/7A 6 = -V D < 0. 
For 7 small, core energy lowering effect can be very large. 

We also notice that the vortex core energy depends on 
a, the distance to the QCP. With the dimensionless quan- 
tity a = a\ 4 /g 2 fi 2 B $l, the change of vortex core energy 

is 5E C ~ —Vo J r ^ xdx(\n 2 x - a) 2 , where r* = Ae _v/ ° is 
the radius where magnetic condensate vanishes. And we 

have 5E C Voe~ 2 ^(3 + 6^ -I- 4a) (see Fig. H. One 

can thus tune the vortex fugacity by changing the dis- 
tance to the QCP. It would be interesting to see whether 
phase diagrams as shown in Fig. [4] can be observed ex- 
perimentally. 

Effect of the magnetic field: In the presence of a per- 
pendicular magnetic field (H _L ab), there will be an 
imbalance of vortices parallel to the magnetic field and 
those anti-parallel, with \n + — n_| > 47J. The un- 
bounded vortices will give rise to finite resistance. When 
the magnetic field is applied parallel to the a6-plane, 
there will be no such effects. This explains the enhanced 
resistivity when applying perpendicular magnetic field 
(Fig. 2c in fl ). One can also see that a small paral- 
lel field will not change Tbkt, i-e. dT/dH c2 \\ — near 
Tbkt, while a small perpendicular field will reduce Tbkt, 
i.e. dT/dH c2 ± < near Tbkt, as observed in Fig. 4a 
of pp. Near Tbkt, where both H c2 \\ and H c2 ± approach 
zero, the ratio H c2{l /H c2± = (dT/dH c2± )/(dT/dH c2ll ) 
thus diverges, as seen in Fig. 3b of pQ. 

Conclusions: In conclusion, we have proposed that su- 
perconducting transition in the heavy fcrmion superlat- 
tice of Mizukami et al. pQ is controlled by BKT transition 
of vortex-antivortex (un)binding. We have also shown 
that magnetic fluctuations modify the conventional BKT 
discussion since they reduce the vortex core energy, and 
thus quantum criticality may strongly influence the phase 
diagram of the vortex system. We made suggestions to 
further test our proposal: The most clear signature of 
the BKT transition is a jump in the superfluid density at 
the transition [34], which can be detected by measuring 
the penetration depth. CeCoIns sandwiched with insu- 
lating layers may make an even better two dimensional 
superconductor. In the opposite limit of a very thin nor- 
mal YbCoIn 5 layer, we expect the crossover to conven- 
tional 3D superconducting transition that also would be 
interesting to test. In a dense vortex matter, vortex- 
antivortex pairs may crystallize, and subsequent melting 
may lead to intermediate hexatic phase [43 44 . It would 
be interesting to look for such phases in systems close to 
a magnetic QCP, where vortex core energy can be sub- 
stantially reduced. 

Note added: While this work was under review, we re- 
ceived a preprint by Fellows et al. [48] , where they study 
a related problem of BKT transition in the presence of 
competing orders, focusing on the behavior near the high 
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Supplementary Material 

A. The dielectric constant 

Here we elaborate on the understanding of the dielec- 
tric constant e c . In BKT theory, the vortex system is 
descibed by the Hamiltonian 

= - ttK [ d 2 r [ dVn(rWr') log |r ^ 1 

- logy J d 2 rn 2 (r), (7) 

where the stiffness K = n s h 2 /ArnksT and the vortex 
fugacity y = e~ E "^ kBT obey the renormalization group 
(RG) equations [35150] 

^K-'il) = 4^y 2 (l), 

|y(0 = [2-KK(l)}y(l). (8) 

Here I = ln(r/£) is the RG scale, £ is the coherence 
length, and E c is the vortex core energy. 

One can define a scale-dependent dielectric constant 
e(r) = K(Q)/K(l), which measures the renormaliza- 
tion of the stiffness K due to the screening of vortex- 
antivortex pairs. Without screening, K takes the bulk 
value K(0) = $gd/167r 3 Ag(T)fc B T, with A b the bulk 
penetration depth. Including the effect of screening, K 
changes with the scale r. One of the most important 
experimental consequencies of the BKT theory is that, 
at the BKT transition temperature, the renormalized 
K 7 i.e. K{1 — oo), approaches a universal value [34) . 
which can be read out directly from the above RG equa- 
tions to be K(oo) = 2/tt. At T = Tbkt,t = oo, the 
scale-dependent dielectric constant becomes of the form 
e(r = oo,T BK t) = <I> 2 d/32TT 2 \ 2 (T BKT )k B T BKT = e c . e c 
is a nonuniversal number. It takes different values for dif- 
ferent systems. For conventional superconductors, e.g. 
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FIG. 5: (a): The dielectric constant e c as function of the 
dimensionless vortex core energy C. The dashed line is 
a fit to the power law behavior. (b): Renormalization 
of the dielectric constant e(r), for different temperatures 
T = Tbkt,0.95Tbkt,0.9T B kt (from top to bottom). Here 
e c = 90, C — 0.0599. (c): The ratio of vortex core energy 
and BKT transition temperature as function of the dielectric 
constant, E c /k B T B KT = {A 1/e /27r)e~ (1 ~ 8)/e , with 9 = 0.83. 



InO x , it is typically 1.1 to 1.9. For YBa 2 Cu 3 07 thin 
films, it is much larger, e c ~ 4.6 [5] or 6 [7] . 

The penetration depth is correspondingly renormalized 
with respect to the bulk value, with A~ 2 = A^ 2 /e(r = 
oo). At the transition, the renormalized penetration 
depth satisfies the relation [34] k B T BKT = <fr 2 l d/32n 2 \ 2 
(Eq. (4) in the main text), which is universal in the sense 
that, different from e c , this relation is identical for dif- 
ferent systems. Thus to determine whether a supercon- 
ducting transition is of the BKT type, it is crucial to 
measure the penetration depth A, and to check whether 
such universal relation between A and Tbkt is satisfied. 
Such relation has been observed in superfuid helium thin 
films [5Tj . 

We can parameterize the vortex fugacity in term of a 
dimensionless quantity C, with y(0) — exp[— CK(0)/A] 
[BJ. C is directly proportional to the vortex core en- 
ergy, with E c = E a C and E = ^d/64n 3 X 2 = 
(e c /27r)fcsTBKT- The vortex core energy can be written 
as E c — (Ce c /2n)k B T BKT . From the above RG equa- 
tions, one can see that the renormalized fugacity vanishes 
at the transition, i.e. y(r — oo,Tbkt) = 0. 

Now we proceed to quantify the relation between the 
vortex core energy E c (or its dimensionless counterpart 
C) and the dielectric constant e c . With the initial con- 
dition K(0) = 2e c /ir, y(0) = e- CK< - ^ 4 and the final 
condition K(oo) = 2/tt, y(oo) = 0, we can numerically 
solve the RG equations. We find that e c = 2,4.6,6,90 
corresponds to C = 7.27,2.24,1.583,0.0599 respectively 
(see Fig. ga)). Following the RG flow (Fig. [gb)), one 
can see that, only very close to the transition temper- 



ature, the dielectric constant changes substantially with 
scale. When moving away from Tbkt, quickly settles 
down to its infared value e^, and decreases signifi- 
cantly with decreasing temperature [BJ. 

It is interesting to notice that for e c > 5, e c and C 
has a power law scaling, e c ~ AC~ e , with the coeffi- 
cient A ~ 8.62 and the power 9 ~ 0.83 (see Fig. |5ja)). 
The dielectric constant and the vortex core energy thus 
has the relation e c ~ A(E C / Eq)~ 9 . A large dielec- 
tric constant corresponds to a small vortex core en- 
ergy. For e c = 90, C = 0.0599, the vortex core energy 
E c = (Ce c /2Tr)k B T BKT ~ {2.7/n)k B T BK ^. Taking 
Tbkt - 1-6^, one obtains E c ~ 0.13meV. For YBCO 
thin films [7j, we have E c ~ (1.583 x 6/2?r) x 7meV ~ 
10.6meV, which is one order of magnitude larger than 
that of heavy fermion superlattice [1] . 

For large e c , we have E c /k B T B KT — 
{A 1 ' 6 /2ir)e~ {1 ~ e)/e (see Fig. |[c)). Due to the small 
power (1 — 9)/ 9 ~ 1/5, for a given Tbkt, a small change 
in the vortex core energy leads to significant change in 
the dielectric constant. Increasing e c from 5 to 90, the 
vortex core energy only changes from 1.54fcgTBKT to 

0.85fc B TBKT- 

In the presence of competing orders, the vortex core 
energy is reduced, E c = Ec — \SE C \. As shown in the 
main text, \6E C \ increases as one approaches the QCP. 
The dielectric constant becomes a function of the dis- 
tance to the QCP, 



= A 



E e V-V n e-> 



v^(3- 



6y/a + 4a) 



(9) 



where a = a\ 4 / g 2 Ms'&o an( ^ a ^ s ^ ne distance to the QCP. 
Vo and a depends on the material specific parameters g, 7. 
In order to determine quantitatively the evolution of the 
dielectric constant near the QCP, more material specific 
microscopic calculations are needed. 



B. Effect of the interface 

At the interface, the Yb ions disorder (due to cross dif- 
fusion and displacements) and act as nonmagnetic impu- 
rities to locally suppress superconductivity in CeCoIn 5 
layers [52] . The superconducting order parameter is 
strongly suppressed near the impurity sites, and it re- 
covers the bulk value over the distance on the order of 



the coherence length [53H5S], £(T) ~ v£ / ^/l - T/T c0 , 
with T c o the bulk superconducting transition tempera- 
ture, £0 the BCS coherence length, and v a number of 
order unity. When the thickness of the CeCoIns layers 
is large, d > £(T), the areas of defect-depressed order 
parameter do not overlap, and the gap is not affected by 
the defects. When the thickness of CeCoIns layers be- 
come smaller than £(T), the depressed areas will start 
to overlap, and the superconducting gap in the CeCoIn 5 
layers will be suppressed. 



G 



At low temperatures with T <C T c o, (,(T) is of order £oj 
which is about the thickness of four layers of CeCoIns. 
So we expect that for n ^> 4, gap has the same value 
as the bulk material; while for n < 4, gap gets sup- 
pressed. This explains the experimental observation that 
the Pauli-limited upper critical field, which is a direct 
measure of the gap, retains the bulk value for n — 5, 7, 
and is suppressed for n = 3. 
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FIG. 6: Illustration of the effect of Yb ions as pair breaking 
nonmagnetic impurities for d > £ and d < £. In the shaded 
regions of size the coherence length £ around the Yb ions, 
superconductivity is suppressed. 
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